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ABSTRACT 

We construct a family of models for the evolution of energetic particles in the starburst galaxy M82 
and compare them to observations to test the calorimeter assumption that all cosmic ray energy is 
radiated in the starburst region. Assuming constant cosmic ray acceleration efficiency with Milky 
Way parameters, we calculate the cosmic-ray proton and primary and secondary electron/positron 
populations as a function of energy. Cosmic rays are injected with Galactic energy distributions 
and clcctron-to-proton ratio via type II supernovae at the observed rate of 0.07 yr _1 . From the 
cosmic ray spectra, we predict the radio synchrotron and 7-ray spectra. To more accurately model 
the radio spectrum, we incorporate a multiphase interstellar medium in the starburst region of M82. 
Our model interstellar medium is highly fragmented with compact dense molecular clouds and dense 
photoionized gas, both embedded in a hot, low density medium in overall pressure equilibrium. The 
spectra predicted by this one-zone model are compared to the observed radio and 7-ray spectra of 
M82. x 2 tests are used with radio and 7-ray observations and a range of model predictions to find 
the best-fit parameters. The best-fit model yields constraints on key parameters in the starburst zone 
of M82, including a magnetic field strength of ~250 ^G and a wind advection speed in the range of 
300-700 km s" 1 . We find that M82 is a good electron calorimeter but not an ideal cosmic-ray proton 
calorimeter and discuss the implications of our results for the astrophysics of the far infrared-radio 
correlation in starburst galaxies. 

Subject headings: galaxies: individual (M82), galaxies: starburst, cosmic rays, gamma rays: theory, 
radio continuum: galaxies 



1. INTRODUCTION 

The coupling between star formation activity and in- 
terstellar gas dynamics is fundamental to galactic struc- 
ture and evolution. Magnetic fields and cosmic rays are a 
key part of that feedback. In this paper, wc fit radio and 
7-ray spectra to study cosmic rays and magnetic fields 
in the archetypal starburst galaxy M82. 

With the discovery of the far infrared-radio luminosity 
(FIR-radio) correlation it became clear that the energetic 
cosmic ray particle population in ga laxies is closely re - 
lated to the star formation process dHelou et al.lfl985h . 
Since far infrared emission from galaxies is primarily 
powered by dust absorption of photon radiation from 
young, massive stars and the 20-cm radio continuum 
by relativistic cosmic-ray electrons, the FIR-radio cor- 
relation requires physical connections to exist between 
star formation processes, cosmic ray production, and the 
physical state of the ISM. These connections have now 
been further extended by 7-ray observations of galaxies 
obtained at TeV energies from the ground with VERI- 
TAS and HESS and from space at lower energies with 
Fermi. 

Starburst galaxies are particularly interesting for test- 
ing models of the FIR-radio correlation in galaxies. 
While electron losses in normal galaxies occur via syn- 
chrotron radiation in a generally diffuse ISM filling large 
volumes, the situation in starbursts is much more ex- 



treme. That this correlation applies to starbursts, with 
their intense star formation, high pressure ISM, and 
large-scale galactic winds, as well as normal galaxies, 
is puzzling as the radio synchrotron emission depends 
on the el ectron ene rgy density and the magnetic field 
strength (jVolk 1989|), as well as losses via galactic winds. 

The observed validity of the FIR-radio correlation has 
led to the concept of calorimeter models for cosmic ray 
interactions in galaxies. The underlying principle of the 
calorimeter model is that the energy radiated by rela- 
tivistic particles scales with the relativistic particle en- 
ergy input. As cosmic rays are produced by supernovae, 
the total energy in cosmic rays follows from the super- 
nova rate and thus from the star- formation rate (SFR). 
A galaxy is then a calorimeter if essentially all the power 
that supernovae transfer to relativistic particles can be 
balanced by observable radiative losses. If the calorime- 
ter model applies to cosmic-ray electrons, then the syn- 
chrotron luminosity is proportional to supernova rate and 
thus the SFR, and the FIR-radio correlation is obtained. 
If it applies to cosmic-ray protons, then the 7-ray lu- 
minosity also depends only on the supernova rate. The 
actual situation is complicated by additional factors, in- 
cluding the effect of the magnetic fields on synchrotron 
radiative efficiencies, non-radiative losses such as ioniza- 
tion, and advection of particles and magnetic fields if a 
galactic wind is present. The existence of such secondary 



2 



Yoast-Hull, Everett, Gallagher & Zweibel 



factors that especially influence synchrotron luminosities 
make the existence of the FIR-radio correlation all the 
more remarkable. 

Observations with HESS and VERITAS on the ground 
and Fermi in space are filling in the gap in measure- 
ments of high energy 7-ray photons from galaxies, in- 
cluding nearby starburst systems (e.g..lAcero et al.ll2009t 
lAbramowski et al.|[2012h lAckermann et al.ll2012l ). Anal- 
ysis of these observations provide unique insights into 
interactions between cosmic rays and the interstellar 
medium (ISM), which in turn lead to a better under- 
standing of the distributi on of feedback ener g y sup- 
plied by supernovae (e.g.. IPaglione et all 119961: [Torres! 



TABLE 1 
YEGZ Model Parameters 



2001 Ide Cea del Pozo et all l2009al lLacki et all 120101: 



Paglione k Abrahamsll2012l ). Further progress in assess 



ing the high energy particle content of galaxies comes 
from radio measurements of the synchrotron emission 
produced by relativistic electrons. Radio observations 
also are advancing in terms of dynamic range, improved 
angular resolution, and sensitivity, especially at low fre- 
quencies. Thus a combination of 7-ray and radio spectra 
can be analyzed to yield information on the high en- 
ergy particle content a nd magnetic fields in galaxies (e.g., 
IThompson eTaII[2006l ). 

We selected M82 for this initial study as this galaxy 
has the advantages of a well determined supernova rate, 
an extensively studied starburst region that includes es- 
timates of the mass of interstellar gas, a carefully stud- 
ied galactic wind, and measurements of the 7-ray and 
radio spectral energy distributions. M82 therefore has 
been a target of choice in testing proton calorimeter mod- 
els to determine the degree to which c osmic ray energy 
is deposited within the galaxy (e.g.. IThompson et all 
2001 IPersic et all l200l Ide Cea del Pqzo et all l2009at 
Lacki et a,l.ll2011b IPaglione k Abrahamsll2012( l 

In this paper we present a one-zone model for calcu- 
lating the electromagnetic power radiated by relativistic 
cosmic-ray protons and electrons in starburst galaxies. 
Our approach is based on analytic calculations applied 
to a starburst zone, where, due to the compact nature 
of the region and high space density of supernovae, cos- 
mic ray diffusion is not important. We include physical 
processes associated with cosmic ray interactions with 
interstellar matter and synchrotron radiation produced 
by primary and secondary cosmic ray electrons. This 
model also incorporates a multi-phase ISM in approxi- 
mate pressure equilibrium and allows for the effects of 
galactic winds that advect particles out of the starburst 
zone. Emergent fluxes are calculated for a model of M82 
and fits to both the radio and 7-ray spectra are used to 
determine key parameters that are not otherwise directly 
constrained by observations of M82. Our approach has 
the advantage of computational simplicity while retain- 
ing the key physics. This allows us to efficiently search a 
wide range of parameters as discussed in §2 to find mod- 
els that best fit the observed radio and 7-ray spectra. 

Our approach to modeling the starburst zone builds 
on that in previous models (e.g.. Ide Cea del Pozo et all 
12009a! : IPaglione k Abrahainsll2012*l ) and explores the ef- 
fects of a wind on the system in detail. Unlike previ- 
ous models, we treat the starburst core and the halo 
separately when dealing with radio emission from the 
galaxy as there is an expected turn-over o f the radio flux 
at low frequencies in the starburst core ( Adebahr et all 



Physical Parameters 



Values Adopted Reference 



Radius SB 
Scale Height SB 
Distance 

Molecular Gas Mass 
Ionized Gas Mass 
Ionized Gas Temperature 
Hot Gas Temperature 
Hot Gas Density 51 
Average ISM Density 51 
IR Luminosity 

Radiation Field Energy Density 51 

SN Explosion Rate 

SN Explosion Energy b 

SN Energy Transferred to CR b 

Ratio of Primary Protons 

to Electrons (N p /N e ) 
Slope of Primary CR 

Source Function 



200 pc 
100 pc 
3.9 Mpc 
2-4 x 10 s M Q 
8xl0 6 M 
8000 K 
6xl0 6 K 
0.33 cm -3 
280-550 cm~ 3 
4xl0 10 L 
1000 eV cm -3 
0.07 yr- 1 
10 51 ergs 
10% 
50 

2.1/2.2 



1 
1 

2 

3,4 
5 
5 
6 
6 



References. [1] iForstcr-Schrcibc r et al.l H2003al): 

Sakai & Mado rH tl99Sfi; [31 INavlor et al.l 120101 ); 14] | Wild et 
I1992I) - | 5] [Fo rstcr-Schrcibcr et al.l (120011); 1 61 IStevens et 
120*03) : [7] IRice et al.l l|19881 ); [8j IFenech et all [|2008h ; 

a Derived from above parameters 
b Excludes neutrino energy 



12013. in pressl ) . Additionally, we use a fully consistent 
energy loss lifetime to test the calorimeter model while 
previous explorations of the cal orimeter model are lim - 
ited to semi-empirical formulae (jAckermann et aDl2012l) . 

The next section describes the properties of M82. Sec- 
tion 3 shows how the population of energetic particles 
was computed. Section 4 contains the details for how 
the model was applied and describes our results and the 
landscape of parameter space. In Section 5, we discuss 
the calorimeter model and implications for the FIR-radio 
correlation, and in Section 6, we present concluding re- 
marks. 



2. PROPERTIES OF M82 
2.1. Structure 



M8 2 is a nearby (D=3.9 Mpc; iSakai k Madord 
119991 ) and well-studied example of a central star- 
burst occurring in a moder ate mass disk galaxy (e.g., 
lO'Connell k Manganolfr978l) . The galaxy is nearly edge 
on, so that much of the starburst zone is optically ob- 
scured by the dusty ISM within the disk. This orien- 
tation, however, allows the galactic wind to be clearly 
seen in wavelengths extending from the x-rays to the ra- 
dio. M82 is a member of the M81 galaxy group and is 
connected to M81 and its surroundings by t idal debris 
that i s readily visible in H 1 2 1-cm line maps (|Yun et all 
119941: iChvnoweth et al]|2008l) . This structure is a result 
of a close passage about 0.2 Gyr ago between M82 and 
M81 that is the probable trigger of the M82 starburst 
(|Yun et all 11991 . The presence of extraplanar H 1 sur- 
roun ding M82 is likel y to also be a result of this interac- 
tion (|Yun et al.lll993h . 

The stellar body of M82 is a high surface bright- 
ness disk with a moderate radial scale length that is 
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only mildly distorted in its outer regions (Ichikawa ct al. 
11995ft . The kinematics of the central region are domi- 
nated by the presence of a bar that extends over much 
of th e ~400 pc diamet er starburst region (jWills et al] 
Il997t iGreve et al.l I2002T ) and is surrounded by a molec- 
ular gas ring or tightly wound arms. This region con- 
tains the majority of the m olecular mass of ~ 3 x 
10 s M 9 (jNavlor et al] 120101 ). which is com parable to 
the m ass of H I surrounding the galaxy (|Yun et all 
119931 ). The star- formation rate within the starburst 
region of ~10 yr" 1 (e.g ., iGao fe Solomon! 120041: 
iForster-Schreiber et al] l2003bD exceeds that of the en- 
tire Milky Way. M82 thus should exhaust its internal 
gas supply in —0.1 Gyr, which is significantly less than 
the time scales of typical galaxies. 

M82's interstellar medium consists of compact, dense 
molecular clouds and dense warm, ionized gas, both em- 
bedded in a hot, low densi ty medium in approximat e 
overall pressure equilibrium (jWestmoquette et al.ll2009f ). 
Pressures in the starburst region are extremely high, on 
the order of 0.5-1.0x 10 7 K cm" 3 (lO'Connell fe Mangand 
[T978tlSmith et al.ll2006l:[Westmoauette et al.H2007l) . Ob- 
servations of CO lines suggest molecular gas densities 
of 10 3 - 10 4 cm" 3 with kinetic tem peratures of 40 K 
([Wild et al.l Il992t iMao et al.l [20001 ) . whereas triatomic 
molecules such as HCO and HNC naturally suggest 
higher densities between 10 4 and 10 48 cm" 3 with tem- 
peratures between 5 and 500 K (e.g..lMuhle et alj|2007t 
iFuente et all 120081: iNavlor et all l2010h~ CO measure- 
ments also constrain the mass of the molecular gas t o 
3 ± 1 x 10 8 M Q (INavlor et aljl20"lfl iWild et al.lfl992h . 
much of which may be in the form of fragmented clouds. 
Observations in the near-infrared give warm, ionized 
gas densities on the order of ~ 10-600 cm" 3 with a gas 
mass of ~ 8 x 10 6 jForster-Schreiber et "all 120011: 
iWestmoauette et al.ll2009h . 

Due its intense star formation, M82 produces its bolo- 
metric luminosity of ~ 4 x 10 10 L© in only a small re- 
gion. Assuming a cylindrical geometry with a radius of 
200 pc and scale height of 100 pc, we find an internal 
radiation field energy density of U ra d ~ 10 3 eV cm" 3 in 
the starburst zone. Thus both the thermal pressure and 
energy density arc substantially enhanced in M82 rela- 
tive to these quantities in the ISM near the Sun. If the 
magnetic field energy density also is increased by factors 
of ~ 10 3 over that in the solar neighborhood, then we 
would expect to find magnetic fields in M82 with aver- 
age strengths of B ~ 150 //G. 

2.2. Galactic Wind 
M82 is well known for its highly visible galactic wind, 



(e.g.. lO'Connell & Mangand 11978 


: iSeaauist & Odeeardl 


19911: Ohvama et al.l 


2002: 


Stevens et al.1 12003: 


Engelbracht et al.l 120061; 


Strickland & Heckmanl I2007D. 



At most wavelengths the outflow has an approximately 
bipolar structure that emerges nearly perpendicular to 
M82's stellar disk. Optical observations of emission lines 
from ionized gas entrained in the wind indicate outflow 
speed s of -500-600 km s" 1 (|Shopbell fc Bland-HawthoTnl 
119981 and references therein). Fits to the x-rays, how- 
ever, indicate higher t erminal wind speeds o f up t o 
-1400-2200 km s" 1 (jStrickland fc Heckmanl [2009h . 
Traditionally the M82 wind has been interpreted as a 



hot outflow powered by s hock heating from supernovae 
(|Chevalier fc Cleeel [19851) , but cos mic ray also could 
be a, factor in driving the wind fereitsch werdt et al] 
119931 ). The variation in speeds and excitation levels 
within the M82 wind indicate that it has a complex 
structure, and so may also be subject to a variety of 
driving mechanisms. 

2.3. Supernova Rate and Cosmic Ray Injection 

An advantage of M82 for cosmic ray studies is that the 
supernova rate can be directly determined through radio 
studies of the numbers and expansion rates of supernova 
remnants. Although values as high as Vsn = 0.2-0.3 
SN yr" 1 have been used for t he supernova rate in M82 
(|de Cea del Pozo et al.ll200"9bl ). we adopt t he more con- 
servat ive estimate of vsn = 0.07 yr" 1 from lFenech et ahl 
(|2008| ). This is likely to be a lower bound on the true 
rate, and is consiste nt with other estim ates of the M82 
star formation rate (|Fenech et al.|[2010l ). 

Supernovae are the assumed accelerators of cosmic 
rays. Of the energy resulting from the explosion, only a 
small fraction is transferred to cosmic rays. In the Milky 
Way, the typical value for t his efficiency factor, 77, is 10% 
(jBlandford fc Eichlerlll987l ). We discuss models for 77 = 
0.04-0.2 in §5. The majority of the energy transferred to 
cosmic rays goes into protons and only a small fraction 
goes into electrons. The ratio of p rotons to electrons is 
here assumed to be N p /N e — 50 (jBlandford fc Eichlerl 
119871 ). The spectrum of cosmic rays accelerated by su- 
pernovae is generally assumed to be a power law with 
spectral index p — 2.1 — 2.2. 

2.4. Radio Continuum 

Of the available ra dio data sets , a col lection of single- 
dish observations by Klein et al.l ([1988D has the largest 
range of frequencies (from 10 7 to 10 14 Hz). Thermal dust 
emission becomes important starting at — 10 11 Hz. The 
lowest frequency data will likely overestimate the flux of 
the starburst due to a large beam size. Because of this 
we only consider data from 1 s to 10 11 Hz in the Klein 
et al. data set. More recently, IWilliams fc Bowerl ([2010D 
observed the starburst region from 1 to 7 GHz with the 
Allen Telescope Array. These interferometer measure- 
ments have significantly smaller fluxes than the single- 
dish observations: the interferometer fluxes are ~12% 
smaller (see Section 4). 

2.5. Gamma Rays 

Starburst galaxies were anticipated to be 7-ray sources 
because of their high star-formation and supernova 
rates. Recently, 7-rays above 700 GeV were detected 
from M82 with the ground-based Cherenkov telescope 
VERITAS. Measurements were made between ~0.9 and 
~5 TeV wi th a flux upper limit (99% confidence level) at 
-6.6 TeV (jAcciari et al.l 120091 ) . M82 has also been de- 
tected in 7-rays with the Fermi "/-Ra y Space Telescope at 
energies from 200 MeV to 300 GeV (jAbdo et al.ll2010[ ). 

3. THEORETICAL APPROACH 

To construct a model that predicts both the radio and 
7-ray spectra, we must first calculate the cosmic ray spec- 
tra for M82. We begin by assuming a source function for 
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Fig. 1. — Lifetimes for cosmic-ray protons (left) and electrons (right). Dotted lines represent the non-interactive loss (advection) time for 
a wind speed of v a ^ v = 500 km s _1 . Dashed lines represent the energy loss time. Solid lines represent the combined particle lifetime. The 
average ISM density is fixed at n = 550 cm -3 corresponding to a molecular gas mass of 4 X 10 s Mq (see Table 1 for additional parameters). 
The discontinuity in the proton lifetime is due to the sudden turn on of pion production as an energy loss at 1.22 GeV. 



cosmic rays accelerated by supernovae. We use a steady- 
state model and assume that diffusion will be relatively 
unimportant in the starburst regiorQ- Then, we calculate 
the primary cosmic ray spectra by including a variety of 
energy losses and an advective galactic wind. Because 
of the very dense nature of the interstellar medium in 
the starburst region and the relatively uniform distribu- 
tion of supernovae in M82, we expect a high interaction 
rate for cosmic-ray protons. As the dominant energy loss 
mechanism is pion production, we calculate the spectra 
for the secondary pions and their associated decay prod- 
ucts (electrons, positrons, and 7-rays). From these pri- 
mary and secondary cosmic rays, we then predict the 
radio and 7-ray spectra. 

3.1. Primary Cosmic Rays 

The co smic ray trans port equation in a homogeneous 
medium {Longair 2011), with diffusion omitted and ad- 
vective losses added, is 



dN(E,t) 
Ft 



d_ 

dE 



dE 
~dt 



N(E,t) 



Q(E,t) 



N(E,t) 



where N(E, t)dE is the number density of particles with 
energy between E and E+dE at time t, dE/dt < is the 
rate at which a particle's energy changes due to radiation 
and collisions, Q(E, t)dE is the rate at which cosmic rays 
with energies between E and E+dE are injected per unit 
volume, and T a dv, which we assume is independent of E, 
is the rate at which particles are advected out of the 
region. 

We assume Q is independent of t and seek time in- 
dependent solutions of eqn. (TTJ. Although the time 
independent version of eqn. (fTJ) is a linear ODE, with 
solutions that can be written down in closed form, the 
solutions are not very useful for our parameter study be- 
cause the complicated form of dE/dt makes them difficult 
to evaluate. Therefore, we make the approximation 



N(E) « Q(E) t(E), 



(2) 



Cosmic ray diffusion times in the Milky Way are about two or- 
ders of magnitude longer than the advection and loss times typical 
of our models for M82. 



where 



-(E)- 1 ee 



adv 



loss 



is the total energy loss rate and 



E 



Tie, 



dE/dt 



(3) 



(4) 



is the energy loss rate due to radiative and collisional 
processes. Equation © becomes exact for T a( i v /Ti oss <C 
1, and is accurate to better than 10% - 20% in cases 
where Q is a power law in E with index near 2, as is 
thought to be the case for acceleration of cosmic rays by 
strong shocks, the leading theory of cosmic ray origin. 

As supernovae are the assumed drivers of cosmic ray 
acceleration, the source function for cosmic rays must be 
related to the total energy input from supernovae. 



Q(E)EdE 



V^snE 51 
V 



(5) 



where h*sN is the supernova rate, V is the volume of the 
starburst region, 77 is the fraction of the supernova energy 
transferred to cosmic rays, and E$i — 1 is 10 51 ergs, 
the typical energy from a supernova explosion. Then, 
assuming the source function is a power law of the form 
Q(E) oc E~ p , the particle spectrum becomes 



N(E) = 



(p - 2) i]iy S NE 5i 



E, 



-p+2 



V 



E~ 



-(E). 



(6) 



Energy losses for protons include ionization, Coulomb in- 
teractions, and pion production. For electrons, losses in- 
clude ionization, bremsstrahlung, inverse Compton scat- 
tering, and synchrotron radiation. Cosmic ray lifetimes 
for values representative of our best fit models are plotted 
in Figure 1, and energy loss rates are plotted in Figure 
2. Energy loss rates are explained in the Appendix. 

Only a small fraction of the supernova power goes into 
cosmic-ray electrons, as compared to cosmic-ray protons. 
We assume the sou rce functions for protons and electrons 
can be related by ()Torresll200l ) 



N 

Qe(E) - jj-Qp(E). 



(7) 
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Fig. 2. — Rate of energy loss for protons (left) and electrons (right). Proton losses include ionization (dashed blue line), the Coulomb 
effect (dotted purple line), and pion production (solid red line). Electron losses include ionization (dot-dashed green line), bremsstrahlung 
(dotted blue line), inverse Compton scattering (dashed purple line), and synchrotron emission (solid red line). A magnetic field strength of 
B = 275 fiG was selected for this example (see §4.4) with the fixed average density of n = 550 cm -3 and a radiation field energy density 
of U rad = 1000 eV cm- 3 . 



We assume the ratio of protons to electrons to be 
N p /N e = 50. 

3.2. Secondaries 

In this section, we discuss pions and the particles re- 
sulting from their decay. Above ~ 1 GeV, the dominant 
energy loss for cosmic-ray protons is pion production, 
and so in this regime, the source function for pions can be 
calculated from the proton spectrum. Assuming a delta 
function approximation, which is appropriate for smooth 
power-law distributions o f cosmic rays, the s ource func- 
tion for pions is given by (|Kelner et al.l 120061 ) 



tion as 



q n (E„) = cn / S(E V - K^T p )a pp {E p )N p {E p )dE p . (8) 



Integrating over proton energy, this becomes 



cn 
n—a p 

Kir 



m p c + 



E-n 



N v m v c + 



E-„ 



(9) 

where K n ~ 0.17 is the mean fraction of the kinetic en- 
ergy of the proton, T p , transferred to the pion per colli- 
sion. 

While the cross sections for pion production from 
proton-proton interactions have been widely studied, 
there is no clear agreement as to which parameteriza- 
tion should b e used. Of th o se ava ilable, we considered 
the model of iKamae et al.l (|2005l ) which includes the 
effects of the diffraction dissociati on process a n d sca l- 
ing violations, and the model of iKelner et al.l (|2006f ). 
whose energy weighted total cross section is based on 
the numerical simulations of proton-proton interactions 
by the SIBYLL code. Comparing the two cross sec- 
tions, the models disagree at energies between 1 and 
10 GeV. However, the 7-ray source functions calculated 
from the two models are nearly identical. The differences 
in cros s section add less than a factor of ~2 in uncer- 
taintv (|Domingo-Santamaria &: Torre s 2005; Lacki et al 

mAs the IKelner et al.l (|2006l ) and IKamae et al 
cross section s give such similar results, we chose 
t o use | Kelner et alls model for simplicity. 
IKelner et al.l (|2006l ) give the total inelastic cross sec- 



<r me i(E p ) = (34.3 + 1.88L + 0.25L 2 ) 

2 

mb, 



(10) 



where L = ln(E p /TeV). This gives us the cross section 
for neutral pion production. For charged pions, we can 
multiply the inelastic cross section by the multiplicity 
of the charged pions. To obtain the multiplicity, we take 
the ratio of the inclusive cross sections for the production 
of charged pions to the inclusive cross section for the 
production of neutr al pions. We u se the inclusive cross 
sections modeled in iDermerl (|1986f ) for this purpose. 

Charged pions are relatively short lived and quickly 
decay into muons and then into electrons, positrons, and 
neutrinos, tt + — > fi + + and /i + — > e + + u e + V p . Be- 
cause the muon moves with nearly the same speed as 
the pion, in the laboratory frame, their source functions 
are equivalent, 9^,(7^) — 9^(7^)- Then, the secondary 
electron and positron source functions are given by 



9e(7e) 



7 ° ( """V; Pile) 



din 



7m 



where 7^ = 7^ ± v^e = ^\fW 
and the electron/positron distribution in the muon's rest 



(11) 

1, y e (max) = 104 



frame is (see lSchlickeiserll2002[) 

27? 



7g 3 (max) 



3- 



27^ 



7g(max) 



Primary and secondary electron and positron spectra can 
be seen in Figure 3. 

Neutral pions have extremely short lifetimes and 
quickly decay into 7-rays, ir° — > 7 + 7. The source func- 
tion for 7-rays is given by (lSchlickeiserll2002l: lPohllll994D 



q 1 ,T T (E. 



1 1 



< Emin VW^ 

where E min = £ 7 + (m 7r c 2 ) 2 /(4£' 7 ). 



--dE r , 



(12) 
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Fig. 3. — Left: Total lifetimes for cosmic-ray protons and electrons. Right: Spectra for primary cosmic-ray protons and electrons and 
secondary cosmic-ray electrons and positrons. Dot-dashed red lines represent protons, solid blue lines represent primary electrons, dashed 
green lines represent secondary electrons, dotted purple lines represent secondary positrons. Parameters were set at B = 275 fiG, v a d v = 
500 km s" 1 , with n = 550 cm" 3 and p = 2.1. 



3.3. 7 -Rays 

In addition to production by the decay of neutral pi- 
ons, 7-rays are also produced by inverse Compton scat- 
tering and by relativistic, non-thermal bremsstrahlung 
from cosmic-ray electrons. In the case of production 
by bremsst rahlung, the source function for the 7-rays 
is given by (lSteckerlll97ll) 



?<y,&rem(-Ky) = CTl(J brem E I N e (E e )dE e , (13) 

JE 1 

where abrem — 3.38 x 10~ 26 cm 2 and N e {E e ) represents 
the combined cosmic-ray electron/positron spectrum. 



For inverse Co mpton scattering, the 7-ray source func- 
(jRvbicki fe Lightmadll979l ) 

Jo e J 7min 7 

(14) 



QyjeiEj) 



16tt J a 

with (|Schlickeiserll2002t) 



In 



(2m e c 2 ) 



1 



1 + 



2 4 
m e c 

~7eZ 



1/2- 



where E 1 is the energy of the resulting 7-ray, e is the 
energy of the incident photon, 7 is the energy of the elec- 
tron. Here, the electron spectrum, n e (^f), is in units of 
cm" 3 . The function F(q,T) is part of the Klein-Nishina 
cross section and is given by (iBlumenthal fc Gouldfl97(ih 



F(q,T) = 2q\n(q) + (l + q-2q 2 ) 
where 



rV(i-g) 
2(1 + r q ) ■ 



r 



4e7 
(m e c 2 ) 



and q 



E~,, 



T("fm e C 2 — Ery 



For the blackbody spectrum, v(e), we use an isotropic, di- 
luted, modified blackbody spectrum (|Persic et al.l 120081 
and references therein) 



where Cda is a spatial dilution factor (given by the 
normalization U ra d = J v(e)ede) and eo corresponds to 
v = 2x 10 12 Hz. 

3.4. Radio Spectrum 

We expect that the majority of the radio flux comes 
from non-thermal synchrotron emission and thermal free- 
free emission. To calculate the radio synchrotr on spec- 
trum, we utilize the emission coefficient given bv lLongaid 

poll . 

dE\ dE 



■synch 



dt 



N(E) 



dv 



(16) 



where N(E) is the combined electron and positron spec- 
trum and 



dE 
~dt 



-a T c 



E 



B_ 2 

2flQ 



T~ ( 17 ) 



gives the synchrotron losses as a function of energy. If we 
make the simplifying assumption that an electron of en- 
ergy E radiates away its energy at the critical frequency, 
v c , such that v ~ v c 



7" • eB/(2irm e ), then 



■synch 
Jv 



-O'T'lTT'eC 



1/2 



2mo 



N(E). (18) 



Because of the warm, ionized gas in the ISM, we ex- 
pect some portion of the radio spectrum to come from 
thermal bremsstrahlung or free-free emission and ab- 
sorption. The free-free em ission coefficient is given by 
(|Rvbicki fc Lightmanl[l97l 



31 



S J = 6.8 x 10" 38 T- 



r 9 ff (19) 

with units of ergs s _1 cm" 3 Hz" 1 . Additionally, we must 
consider the inverse process, free- free absorption. The 
absor ption coefficient is given by dRvbicki fc Lightmanl 
H979h 



Jf _ 



= 0.018T 



-3/2 



Z n e riiU 2 gff cm . 



(20) 

In the "s mall-angle, classical regi on," the mean Gaunt 
factor is (jNovikov fc Thorndll973l ) 
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Fig. 4. — Schematic of the process. The red explosion is a supernova remnant, the black cloud is a molecular cloud, the green line is a 
schematic magnetic ficldlinc loaded with cosmic rays, and the blue arrow represents a galactic wind flow. The figure on the left shows an 
initial configuration, and the figure on the right shows the configuration after the fieldline has been carried out some distance by the wind. 



where Ry = 13.6 eV is the Rydberg constant and ( = 
e 7 w 1.781 with 7 being Euler's constant. 

For both emi ssion and absorption, the ra diative trans- 
fer equation is (jRvbicki fc Lightmanlll979t ) 



(22) 



where S v = j v /K v is the source function. We discuss 
solutions to cqn. ([22]) in §4.3.1. 

4. MODEL AND RESULTS 

We have calculated the steady-state primary and sec- 
ondary cosmic ray spectra according to the assumptions 
made in Section 3. From there, we compared the ob- 
served properties of M82 with our model. We included 
the multiphase nature of the interstellar medium in or- 
der to better predict absorption and emission in the radio 
spectrum. We calculated the synchrotron emission from 
the primary and secondary electrons and positrons and 
we also calculated the effects of free-free emission and 
absorption on the radio spectrum. In this section, for 
both the radio and 7-ray spectra, we discuss the effects 
and constraints of various parameters on the predicted 
spectra and the resulting best-fit models. 

4.1. Cosmic Ray Propagation in the Disk 

We assume that the cosmic rays sample the mean den- 
sity in the disk. This assumption is not trivial, since we 
expect that most of the cosmic rays are accelerated in the 
hot, low density medium, which is flowing out of the disk, 
while most of the mass is in the form of dense, molecu- 
lar gas which has a small filling factor. Cosmic rays will 
sample the mean density if (1) most of the magnetic field 
lines they are accelerated on connect to a "fair sample" 
of the interstellar gas, and (2) cosmic rays can propa- 
gate a sufficient distance along the fieldlines, or diffuse 



sufficiently across the fieldlines, to encounter dense gas 
before they - and the fieldlines - are borne out of the disk 
by the wind (see Figure 4). The timcscale for the latter 
is T a dv, the escape time introduced in eqn. (5). 

By a "fair sample," we mean that the ratio of inter- 
stellar gas mass to magnetic flux is uniform within the 
starburst region. Since the covering factor of the molecu- 
lar clouds is below unity, this suggests that the fieldlines 
have a substantial random component, as is the case in 
the Milky Way. In the absence of any information to the 
contrary, we assume the fieldlines arc indeed uniformly 
loaded, so that the average gas density on any magnetic 
flux tube within the disk is the mean interstellar density, 

fmol^mol fion^ion ~t~ (1 fmol fion)^hot ~ fmol^mol- If 

there are N mc molecular clouds in the disk and the mean 
dimension of the disk is L ~ y 1 ! 3 where V = 2-kR 2 H. 



then the mean distance L to a cloud is of order LN, 



-1/3 



If the supernova remnants where the cosmic rays are ac- 
celerated and the molecular clouds are spatially uncorre- 
lated, then l c is the mean minimum distance that cosmic 
rays must travel to encounter dense gas. 

In order to estimate the time it takes cosmic ray parti- 
cles to travel a di stance l c . we as sume the cosmic rays are 
self-confined (see lKulsru d (2005) for a pedagogical treat- 
ment): any cosmic ray density gradient produced by a 
point source of cosmic rays, or by the global distribu- 
tion of cosmic ray sources within the starburst region, is 
associated with an anisotropy of the cosmic rays in veloc- 
ity space. If the mean cosmic ray velocity, or streaming 
speed, exceeds the Alfven speed va, then the cosmic rays 
destabilize short wavelength Alfven waves. The waves 
grow large enough to scatter the cosmic rays, reducing 
their drift speed to va- In this picture, the propaga- 
tion time is I c /va- The cosmic rays will encounter dense 
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Fig. 5. — 7-ray spectra. Left: 7-ray spectrum with parameters p 
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100 cm -3 . This spectrum is the best fit to the 7-ray data. Right: 7-ray spectrum with parameters p = 2.2. M mo i = 4 X 10 s Mq, B = 
275 AtG, v adv = 400 km s" 1 , 

^ion — 150 cm . This is the 7-ray spectrum that corresponds to the best radio fit for a spectral index of 
p = 2.2. The solid black lines represent the total 7-ray flux, the dashed red lines represen t the contribution from neutral pion decay, and 
the dotted blue line s represent the contribution from bremsstrahlung. 7-ray data include: Ackcrmann ct al. (2012) [Fermi - blue circles), 
lAcciari ct al. (2003) (VERITAS - red squares). Data with downward arrows represent upper limits for both Fermi and VERITAS data. 



clouds if 



VATadv 

Tc 



Vadv 



2/3 



N, 



1/3 



(2tt) 1 /3 



> 1, 



(23) 



where in the last step we have identified r a dv with an 
advection time H/v a( i v . 

In §4.4 we will show that both the radio and 7-ray 
spectra are very well fit by a simple family of models 
in which the magnetic field strength B and advection 
velocity v a dv are related by the empirical formula 



Vadv 



4.0B, t - 600, 



(24) 



where is the magnetic field in fiG and v ac i v is in 
km s" 1 . Equation dHJ) holds for 150 < B^ < 350. Sub- 
stituting this relation into eqn. (|23[) . using H/R = 0.5 
(see Table 1) and solving for N mc , we find that cosmic 
rays will encounter dense clouds before being advected 
out of the Galaxy if 



N mr > 8tt 2.00 - 



300 



3/2 
l hof 



(25) 



According to eqn. ([2~5j) . for the value of N mc needed to 
satisfy the criterion varies from to 7 as B^ ranges 
from 150 to 350 with n/ lot = 0.33 cm~ 3 (see Table 
1). With the molecular gas mass of the M82 starburst 
region estimated to be about 3 x 10 8 Mq, and giant 
molecular cloud masses in the Milky Way estimated at 
10 4-6 M , this criterion would seem to be easily sat- 
isfied. Moreover, the drift speed actually exceeds va 
if the waves are strongly damped, making eqn. (|23j) 
even easier to satisfy. A more careful assessment based 
on a Monte Carlo simulation suppor ts this conclusion 
(|Boettcher et al.ll2013. in preparation!) . This effect may 
be important for the higher energy cosmic rays. Once the 
cosmic rays reach the weakly ionized clouds, the waves 
which scatter them are strongly damped, and the cos- 
mic rays interact primar ily through physical collisions 
(|Everett fc Zweibelll2l)Tl"l ). 

4.2. -y-Rays 



Though we incorporate the multiphase nature of the 
interstellar medium into our model for the radio spec- 
trum (see §4.3.1), this is not necessary in the case of 
the 7-rays because of our assumption that they sam- 
ple the mean density (§4.1) which is almost entirely 
due to molecular clouds. The ionized gas mass is well- 
constrained by observations to be Mi on = 8 x 10 6 M Q . 
In §4.4 we describe models with molecular masses of 
M mo i = 2 — 4 x 10 8 Mq corresponding to mean densities 
of (n) = 280- 550 cm" 3 . 

Pion decay is the main mechanism for the production 
of 7-rays. When advective losses dominate over energy 
losses, the 7-ray flux is proportional to the gas density. 
However, in the calorimeter limit, where cosmic-ray pro- 
ton energy losses dominate over advective losses, the 7- 
ray flux is essentially independent of density and so the 
number of 7-rays depends only on the number of pro- 
tons. Other parameters that affect the 7-ray spectrum 
include the supernova rate (i^sn = 0.07-0.1 SN yr" 1 ) 
and the distance (d = 3.5-3.9 Mpc), both of which are 
observationally constrained. 

In addition to the supernova rate and distance, wind 
advection speed (v a( i v ) affects the 7-ray spectrum. For 
our cosmic ray spectra, we use a combined lifetime (see 
eqn. Q) which includes energy losses and advection. 
The advection time is given by T a d v = H/v a dv where 
H is the height of the starburst region and v a d v is the 
speed of particles in the wind in the starburst region. 
Of the parameters that affect the 7-ray spectrum, the 
least constrained is the wind speed. For the 7-ray spec- 
trum, the main constraint on the wind speed is the Fermi 
observations. The proton lifetime peaks between 1 to 
10 GeV, and it is this portion of the proton spectrum 
that contributes to the energy range for the 7-ray spec- 
trum observed by Fermi (0.3 to 30 GeV). So, if the wind 
speed becomes too large, the maximum proton lifetime 
decreases and the 7-ray spectrum will no longer agree 
with the Fermi observations. If this is the case, then the 
galaxy is no longer calorimetric. 

Assuming values typical of the Milky Way for parame- 
ters related to cosmic ray acceleration by supernovae, we 
find that the 7-ray flux from neutral pion decay fits the 
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TeV energy 7-ray data and all but the two lowest-energy 
Fermi data points as seen in Figure 5 and discussed in 
§4.4. As the 7-ray flux from neutral pion decay does not 
account for all of the observed flux below ~1 GeV, there 
must be another emission mechanism responsible for the 
total 7-ray flux at these energies. Two possible mecha- 
nisms are inverse Compton emission and bremsstrahlung. 
For 7-ray emission from the inverse Compton effect, we 
found that the 7-ray flux in the starburst region was on 
the order of 2 x 10~ 10 GeV cm -2 s _1 at ~1 GeV assum- 
ing an approximate radiation field energy density of 1000 
cV cm" 3 (sec Table I). While our model currently in- 
cludes inverse Compton emission from only primary elec- 
trons, we found that it to be a significant contribution 
from the total 7-ray flux at both GeV and TeV ener- 
gies. Inclusion of emission from secondary electrons and 
positrons will increase the total 7-ray flux by a factor of 
3 at minimum. We will provide a more accurate model 
for inverse Compton emission in future work. 

We also calculated the 7-ray flux from non-thermal, 
relativistic bremsstrahlung. As bremsstrahlung depends 
on the electron spectrum, the flux is dependent upon 
the magnetic field strength and wind speed. We find 
that bremsstrahlung makes a significant contribution to 
the 7-ray flux below ^10 GeV (see Figure 5). Thus, 
to accurately model the observed 7-ray flux at lower 
energies, we must include both neutral pion decay and 
bremsstrahlung in our work. 

4.3. Radio Spectrum 

In the case of the radio spectrum, we have several 
more parameters that affect our fit than for the 7-ray 
spectrum. These parameters are the ionized gas density 
(ni on ), the magnetic field strength (B), and the wind 
speed (v a dv)- Free- free emission and absorption both 
scale as the square of the gas density. Thus, the lower 
the gas density, the lower the frequency at which the 
radio spectrum turns over and the lower the contribu- 
tion of free-free emission to the radio spectrum at higher 
frequencies. Additionally, synchrotron emission is scaled 
by the square of the magnetic field strength (B), and the 
wind speed (v a dv) affects the cosmic ray "dwell time" in 
the starburst. To quantitatively distinguish between the 
effects of our parameters, we use \ 2 tests to calculate the 
goodncss-of-fit. 

4.3.1. Multiphase Solutions for Radiative Transfer 

The emergent radio spectrum is due to free-free and 
synchrotron emission, modified by free-free absorption 
as discussed in §3.3. Fitting a model to the observed 
radio spectral energy distribution therefore depends on 
knowing the degree of free-free absorption. This is most 
pronounced at lower frequencies and these observations 
are essential for obtaining a physically realistic fit to the 
radio data for M82. We tested several models and found 
that the free-free covering factor was always near unity. 
A simpler model consisting of a cylindrical wall of warm, 
ionized gas surrounding a lower density hot medium pro- 
vides a reasonable representation of the free-free radio 
absorption in M82. 

The multiphase nature of the ISM is incorporated by 
considering the different properties of the warm and hot 
ionized gas. The warm, ionized gas is assumed to be mag- 
netically threaded and produces both synchrotron and 



■synch 
Jv 



jff,ion*j ag we u ag 



free-free emission (jwiM 

free-free absorption (kwim = «v ' ion )- Synchrotron self- 
absorption is found to be negligible in M82, so we need 
only consider free-free absorption. The hot, low den- 
sity gas also produces synchrotron, and due to its lower 
density, negligible amounts of free-free emission and ab- 
sorption { Jhot =jy nch +j f J Mt «#""*). 

Adopting our assumed cylindrical geometry with a 
shell of warm, ionized gas, we take the shell width to 
be 1/2 where I is the radius of the cylinder. Because the 
galaxy is nearly edge-on, we see emission from both the 
front and the back of the cylinder. Emission from the 
front of the cylinder is 



JWIM 

kwim 



^2 _ g — kwimI/2^ 



(26) 



Emission from the rear of the cylinder can be absorbed 
within both the front and back warm ionized gas shells. 
Thus, 



' v,back 



JWIM 
KWIM 



^2 g-KWmf/^j ^ — kwimI/2 (27) 



Then, the radiative intensity emergent from the hot, dif- 
fuse gas is given by 



L,hot = 2jhotrie- KwlMl ' 2 , 



(28) 



where r% is the radius of the central volume of hot, diffuse 
gas (r starburst = Ti + 1/2). Then, the total emergent 
intensity from the starburst region is 



JWIM 

kwim 



(l - e ~ KwlMl ) + 2j hot r i e~ K ' WIMl/2 . (29) 



This model ignores the presence of a radio halo. Obser- 
vations of the radio halo shows it has a steeper spectral 
index than that of the main body and thus can become 
dominant at low frequencies (jSeaquist fc Odegardlll991l : 
lAdebahr et al.ll2013. in press! ) . We therefore constructed 
a model that assumes all of the 74 MHz flux from M82 
arises in an unabsorbed synchrotron halo. This assump- 
tion is consistent with the parameters of our cylindrical 
model, where any 74 MHz emission from the main star- 
burst zone is heavily absorbed. We assume a standard 
-0.7 spectral index for the halo. The radio halo's spectral 
energy distribution then is simply 



h. 



halo 



-0.7 



(30) 



where Io,haio is determined by our lowest frequency radio 
data point. The final expression for intensity is 



/„ = m±L (i 



-KWIM l 



)+2j hot r l e- KW ' Ml ' 2 + I uM io. 

(31) 

Figure 6 shows the radio spectrum predicted by the 
simple cylindrical model, which shows the contributions 
from the model components. Emission from the hot, dif- 
fuse gas, eqn. (f28|) . is represented by the dotted blue 
line; emission from the warm, ionized gas, eqn. (|26|) & 
(|2"T|) . is represented by the dot-dashed green line; and 
emission from the halo, eqn. (|30[) . is represented by the 
dashed red line. Note that the free- free high frequency 
emission is under estimated since we cannot use a sin- 
gle value for the mean electron density (see §4.4). Our 
best-fit model predicts a flux from the starburst region 
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Fig. 6. — Best fit to radio spectrum (x 2 = 22.6) with parameters B = 275 /iG, v a d v = 500 km s 1 , ra; on = 100 cm 3 with p = 
2.1, M mo i = 4 X 10 8 Mq. The solid black line denotes total radio flux, the red dashed line represents the radio emission in the halo, 
the dotted blue line represents the radio emission in the h ot, diffuse gas, an d the dot-dashed g r een lin e represents r adio emission in the 
warm, ionized gas. Radi o data include [Klein et al. (198 81) (blue triangles), IWilliams & Bower (2010) (red circles ), Cohen ct al.l 1)20071 ) 
(green circle), Basu ct al. (2012) (green circle), Carlstrom & Kronbcrg (1991) (green circle), Adcbahr et al. (2013, in press) (cyan circle). 
Triangles represent single-dish data and circles represent interferometer data. Open squares represent the median data set used for the x 2 
tests. 



at 333 MHz in a greement with direct measurements from 
lAdebahr et al.l (|2013. in press! ). 



TABLE 2 
Fitted Model Parameters 



4.3.2. Radio Observations 

Because o f the d ifferi ng methods of obser v ation for the 
IKlein et al.l (|1988D and IWilliams fc Bowed (|2010D . data 
sets it is difficult to say which is more accurate. At most, 
we can say that the single-dish observations serve as an 
upper limit while the interferometer observations provide 
a lower limit for the radio fluxes for the starburst region 
of M82. For comparison with our model, we calculated 
the median bet ween the two data sets and scaled down 
the data from IKlein et all (|1988| ) a s it encompasses a 
larger range of frequencies than the IWilliams fe Bowerl 
dat a. We also i ncreased the error bars to six percent for 
the IKlein et al.l data to include both data sets. 

We do not includ e the s ingle-dish 87.2 G Hz measure- 
mcnt of Uura et ail (|1978| ) , included in the IKlein et all 
data set, as it is u nusually low compared to even the in- 
terfer ometer data ([Klein et al.l 119881 : IWilliams fc Bowerl 
[20101) . Instead , we use the 92 GHz interferometer mea- 
surement from iCarlstrom fc Kronberd (|1991l ) . Addition- 
ally, we i nclude newer low frequency meas urements at 
74 M Hz (jCohen et al.ll2007l ) and 333 MHz (jBasu et al.l 
I2012T) . As these are both interferometer observations, we 
scale them up to use with our median data set. Distin- 
guishing between core and halo emission at these low fre - 
quencies is quite difficult (|Adebahr et al.ll2013. in press! ) 
and so we incorporate a halo emission component into 
our model (see above). 

4.4. Results 
AAA. x 2 Tests 



Physical Parameters 



Tested Range Reference 



Magnetic Field Strength (B) 
Advcction (Wind) Speed (v ac i v ) 
Ionized Gas Density (rii on ) 



50-350 nG 
0-2000 km s" 
50-500 cm" 3 



1.2 
3.1 

5 
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One of the primary goals of developing a one-zone 
model is to be able to find a set of parameters which 
produce an acceptable fit to both the radio and 7-ray 
spectra. To find the best-fit combined solutions, we use 
X 2 tests on both spectra. We begin by using observations 
and previous studies to determine the range of values 
to test for each parameter. Previous studies have cal- 
culated th e equipartition magnetic field strength to be 
-150 /iG (|de Cea del Pozo et alJ[2009al ). Thus, we test 
values from 50 to 350 /xG. For the wind speed, model fits 
to x-ray observations of t he wind give wind speeds o n 
the order of 1500 km s" 1 ([Strickland fc Heckmanll2009t) . 
while optical obser vations give wind speeds on the or - 
der of 500 km s" 1 ()Shopbell fc Bland-Hawthornlll998| ). 
The two measurements do not inherently contradict each 
other as they mesaurc different materials in the wind. 
Thus, we test values from to 2000 km s — 1 . We have 
fixed the total mass of the molecular and ionized gases, 
and in doing so, we have also fixed the proton density 
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Fig. 7. — Contour plots showing changes in \ values for changes in magnetic field strength (B) and advection (wind) speed (v a d v ). Top 



left: Parameters set at p = 
= 150 cm" 3 . Bottom left: 
M mol = 4 x 10 8 M Q , n ion 



2.1, M mol = 2 x 10 s 
Parameters set at p 
= 150 cm" 



M S , 
= 2.1, 



,-3 



L ion — 

M mol 



300 cm 3 . Top right: Parameters set at p = 2.1, M mo l = 3 X 10 8 Mq, n 



4 X 10 s Mr. 



Observations constrain the ionized gas densities to be 
10 to 600 cm" 3 with kinetic temperatures of ^8000 K 
(|Forster-Schreiber et al.|[200TI ). We test values from 50 
to 500 cm -3 (see Table 2). Additionally, we test two dif- 
ferent cosmic ray spectral indices (p = 2.1, 2.2) and three 
different molecular gas masses (M mo i = 2 x 10 s , 3 x 10 s , 
4 x 10 s M ). We took advantage of the computational 
simplicity of our model to run 49,000+ cases. 

In determining the best fit to the radio spectrum, we 
performed \ 2 tests using the available radio observations, 
discussed above. Fits to the eleven radio data points 



100 cm" 



Bottom right: Parameters set at p 



I on 

2.2, 



were used to constrain the ionized gas density (n, on ), the 
advection speed (v a dv), and the magnetic field strength 
(B) (see Figure 7). For the 7-rays, fits to the eight data 
points (excluding upper limits) were used to place addi- 
tional constraints on the advection speed and the mag- 
netic field strength. 

In looking for the joint best-fit model, we found that 
for the radio spectrum there is a degeneracy between 
magnetic field strength and wind speed such that as wind 
speed increases, an increase in magnetic field strength al- 
lows a similiarly good fit (see Figure 7). We explain the 
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Fig. 8. — All fits within 3<r of minimum x 2 for n; on = 100 cm -3 with p = 2.1, M moi = 4 X 10 s Mq. The black line represents the 
best-f it (y_ 2 = 22.6) and the grey lines r epresent all ot her 3a fits. Radio d ata include Klein et al. ( 1988) (b lue tr i angle s) , IWilliams &: Bowerl 
1120101 ') (red circles). ICohen et al. (200711) (green circlet. rBasu et al.l ||20T2T ) ( green circle). ICarlstrom Hi Kronbered 1)19911 ) (green circle) . Open 
squares represent the median data set used for the \ 2 tests. 



reasons for this degeneracy below. For the 7-ray spec- 
trum we also find a valley of solutions (several minima 
as opposed to a single absolute minimum). If we restrict 
the possible 7-ray models to those deemed to have ac- 
ceptable radio spectra, we find that the % 2 tests for the 
7-ray spectrum limits the tested wind speeds to only a 
few. 

The best-fit solution has a \ 2 value of \ 2 = 22.6 for 
the radio spectrum and \ 2 = 9-6 for the 7-ray spectrum. 
It occurs for an ionized gas density of 100 cm -3 , a mag- 
netic field strength of B = 275 /iG, and a wind speed 
of v a dv = 500 km s _1 (Table 3). The models within 3a 
have magnetic field strengths of B = 225-350 /xG, ad- 
vection speeds of v a( i v — 300-700 km s , and ionized 
gas densities of rii on = 50-250 cm -3 . Figure 8 shows 
all radio spectra for models within 3<t of the minimum 
for a single spectral index and molecular gas mass. Of 
the 16,000+ models tested, we found 666 radio models 
within 3<7 while only 89 had both radio and 7-ray models 
within la. 

The best-fit 7-ray model can be seen in Figure 5 and 
the corresponsing radio model can be seen in Figure 6. 
The radio spectrum agrees to within the errors for ten of 
eleven data points, while the 7-ray spectrum agrees with 
six of eight data points. It can be clearly seen in Figure 6 
that the radio model does not agree with the 92 GHz data 
point. The flux at 92 GHz has been previously shown to 
be m ostly due to free-free emission (jWilliams fc Bowerl 
l20T0t) . Because we have limited the ionized gas to a single 
density, we achieve the turnover at low frequencies due 
to free-free absorption but not the required flux at high 
frequencies due to free-free emission. To achieve both the 
observed free-free absorption and emission would require 
another level of detail to the model. Thus, in maintaining 
the simplicity of the model, we are not able to produce 
the required free-free flux at high frequencies. 



TABLE 3 
Best-Fit Model Parameters 



Physical Parameters 



Best-Fit Value 



Magnetic Field Strength (B) 


275 




Advection (Wind) Speed (v at i v ) 


500 


km s _1 


Ionized Gas Density (rii on ) 


100 


cm 


Spectral Index (p) 


2.1 




Molecular Gas Mas (M mo ;) 


4 X 


10 8 M Q 



Note. 



Results for y 2 .. = 22.6, y 2 



9.6 



4.4.2. Wind Speed and Magnetic Field Strength 

In the course of our % 2 tests, we found a degeneracy 
between the magnetic field strength and the wind speed. 
As one increases the magnetic field strength, you must 
also increase the wind speed to achieve a similarly good 
fit. Why might such a degeneracy occur? 

The magnetic field strength and the wind speed both 
affect the amount of synchrotron emission. The mag- 
netic field essentially scales the radio spectrum up and 
down as the square of the magnetic field strength; how- 
ever, because the high frequency portion of the spectrum 
is dominated by thermal emission from the ionized gas, 
changing the magnetic field does not change the high fre- 
quency end. If we increase the wind speed, then we de- 
crease the advection time and thus decrease the number 
of cosmic-ray electrons and positrons in the starburst re- 
gion. By advecting away some of the electrons/positrons, 
we decrease the amount of synchrotron emission. Be- 
cause we have fixed the molecular and ionized gas masses 
and therefore the average density, it is necessary to have 
advection to create a better fit as otherwise we ovcrprc- 
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Fig. 9. — Contour plots showing changes in x 2 values for changes in magnetic field strength (B) and advection (wind) speed (v a d v ) 
for acceleration efficiencies of 4% and 20%. Left: Parameters set at p = 2.1, M mo ; = 4 X 10 s Mq, rii on = 150 cm~ 3 , r\ = 0.04. Right: 
Parameters set at p = 2.1, M moi = 4 X 10 s Mq, n; on = 50 cm" 3 , n = 0.2. 



diet the radio spectrum. If we increase the wind speed 
too much, then the advection time is shorter than the 
loss lifetime, and we no longer have enough cosmic-ray 
protons to produce the correct amount of 7-ray flux. 

4.4.3. Gas Mass and Spectral Index 

In addition to varying the values for ionized gas den- 
sity, magnetic field strength, and wind speed, we also 
tested different values for spectral index and molecular 
gas mass. We test a much smaller range of values for 
the cosmic ray injection spectral index (p = 2.1, 2.2) and 
molecular gas mass (M mol = 2 x 10 8 , 3 x 10 8 , 4 x 10 8 M ) 
than for the other parameters. While the observations 
relating to ionized gas density, magnetic field strength, 
and wind speed leave much uncertainty as to their true 
values, observations for spectral index and molecular gas 
mass are more constraining. Thus, we only need test a 
few values for each variable. 

We find that for both the radio and 7-ray spectra, the 
spectral index of p = 2.1 is a better fit than p = 2.2. 
While selecting p = 2.1 over p = 2.2 leads to less than a 
factor of 2 reduction in the minimum \ 2 value, the 7-ray 
models withp = 2.2 consistently miss the TeV VERITAS 
data points (see Figure 5). For molecular gas mass, we 
find that the minimum \ 2 value decreases as the gas mass 
increases with there being a slightly larger gap between 
in x 2 values between 2 x 10 8 and 3 x 10 8 than between 
3 x 10 8 and 4 x 10 8 . 

5. DISCUSSION 

Over the course of this paper, we have investigated a 
model for cosmic ray interactions in the starburst galaxy 
M82. We incorporate synchrotron emission and free- free 
absorption and emission into the radio spectrum and in- 
clude neutral pion decay and bremsstrahlung as mecha- 
nisms for 7-ray production. For the most part, the fits 



agree well with the data, with the exception of 92 GHz 
radio data point. The radio fits also indicate that there is 
a degeneracy between magnetic field strength and wind 
speed. We find that for both the radio and 7-ray spectra 
a spectral index of p = 2.1 and a molecular gas mass of 
M mo i = 4 x 10 8 M produce the best fit (see §4.4) , but 
a variety of different models fit both the radio and 7-ray 
spectra well. However, while the \ 2 tests on the radio 
spectrum allow for fits over a large range of wind speeds, 
the x 2 tests on the 7-ray spectrum more severely limit 
the range of allowable wind speeds. 

5.1. Is M82 a Calorimeter? 

One of the main questions we seek to answer in this 
paper is whether or not M82 is a calorimeter. We define 
a calorimeter such that cosmic-ray proton or electron en- 
ergy losses dominate over advective or diffusive losses. In 
the case of the cosmic-ray electrons, it is clear that the 
galaxy is an electron calorimeter as energy losses domi- 
nate over advective losses for all wind speeds tested (see 
Figure 1). 

However, it is more difficult to decide exactly when a 
galaxy is no longer a proton calorimeter. In the case of 
the cosmic-ray protons, the advection timescale is about 
equal to the peak of the energy loss lifetime (see Fig- 
ure 1) meaning that the galaxy is a 50% calorimeter. 
How calorimetric a galaxy is depends on not only on the 
wind speed, which determines the advection timescale, 
and the average density, which determines the energy 
loss timescale, but also on the cosmic ray acceleration 
efficiency. Throughout the paper, we have assumed cos- 
mic ray acceleration efficiency (r/) of 10%. To investigate 
the impact of acceleration efficiency, we tested two al- 
ternative values: a decreased efficiency of r\ = 4% and 
an increased efficiency of 77 = 20% (that could also be 
achieved with an efficiency of 10% and an supernova rate 
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of 0.14 yr" 1 ). 

For the 4% case, the resulting radio models consis- 
tently required higher magnetic field strengths due to 
a lack of electrons because of the decrease in the pri- 
mary electron population and because of a decrease in 
the secondary electron/positron population because of 
the decrease in primary protons. The best models re- 
quired wind speeds on the order of 0-100 km s -1 (see 
Figure 9). Wind speeds of this magnitude give advection 
timcscales that are longer than the proton energy loss 
timescales, resulting in an 80%-100% calorimeter. The 
best fit had a % 2 value nearly a factor of 2 larger than the 
best fit for an efficiency of 77 = 10% for both the 7-ray 
and radio spectra. For all three acceleration efficiencies, 
we found 317 of the 49,000+ models to be within la for 
both the radio and 7-ray spectra. Of those 317 models, 
only 45 had an acceleration efficiency of 4%. Based on 
these results and the extremely low wind speeds, we can 
effectively rule out the 4% case. 

Conversely, for the 20% efficiency case, we found 
that significantly higher wind speeds were required to 
achieve good fits because of an increase in the over- 
all electron/positron population. The best models re- 
quired wind speeds on the order of ~1500 km s _1 
(see Figure 9). While models for an efficiency of 10% 
agree with measurements of wind speeds in the optical 
(jShopbell fc Bland- Hawtho"rnlll998ft . wind speeds on the 
order of 1500 km s -1 agree w ith x-ray measurements by 
IStrickland fc Heckmanl (l2009l ). 

Comparing the resulting advection timescales to the 
proton energy loss timescales gives a 35% calorimeter. 
We find that for this higher acceleration efficiency, the 
minimum \ 2 value for the radio spectrum is smaller than 
in the 10% case, while the \ 2 value for the 7-ray spectrum 
is just slightly higher than before (x 2 = 10.8 vs \ 2 = 9.6). 
Thus, while we rule out the 4% acceleration efficiency, we 
consider the 20% case to be plausible. 

Acceleration efficiency directly affects our initial cos- 
mic ray population. However, it is not the only factor. 
For example, we use the lower bound of 0.07 SN yr -1 
for our supernova rate. If we scaled the supernova rate 
up and the acceleration efficiency down in equal measure, 
our initial population of cosmic rays would be unchanged. 
While we have selected one particular set of parameters 
for our initial cosmic ray population, it is simple to ob- 
tain the same population from a slightly different set of 
initial parameters and as such, we would again produce a 
degenerate relationship between magnetic field strength 
and wind speed. Thus, we would produce the same fits 
as presented above. 

5.2. FIR-Radio Correlation 

It remains a curious question as to why starburst galax- 
ies lie on the same FIR-radio correlation as the normal, 
quiescent galaxies. The general explanation for the FIR- 
radio correlati on in d isk ga laxies is the electron calorime- 
ter theory by IVolkl (p89ft . This paper postulates that 
both FIR and radio synchrotron emission arc propor- 
tional to the supernova rate as FIR emission is essentially 
starlight reradiated by dust grains and radio synchrotron 
emission co mes from en ergetic electrons accelerated by 
supcrnovae ([Volkl 11989ft . As such, it is required that 
cosmic-ray electrons radiate away their energy wihtin the 



galaxy to achieve the observed radio synchrotron emis- 
sion. However, in dense environments such as starburst 
galaxies, there is a second source of electrons/positrons 
from proton-proton collisions. 

The overall importance of secondary elec- 
trons/positrons is dependent upon the ratio of primary 
protons to primary electrons (N p /N e ). The significance 
of secondary electron/positrons decreases as the ratio of 
primary protons to electrons decreases. In our model, we 
leave this ratio as an independent free parameter fixed 
at 50 based upon supernovae studies in the Milky Way. 
ISchlickeiserl ((2002ft makes a case for this parameter being 
model dependent such that the ratio depends on the 
masses of protons and electrons and the selected spectral 
index, N p /N e = {m p jm^) { -' p ~^l 2 . Our prefered spectral 
index is p = 2.1, which gives a ratio of N p /N e ~ 60. As 
our free parameter selection is smaller than this value, 
the importance of secondary electrons/positrons would 
only increase by making the ratio model dependent. 

Thus, while we find that M82 is an electron calorime- 
ter, the overall picture is complex. It is true that all 
cosmic ray electrons are mainly confined to the starburst 
core. However, because the galaxy is not also a proton 
calorimeter, there is an inherent loss of any secondary 
electrons and positrons that would have resulted from 
proton-proton collisions in the core. Thus, any electrons 
produced in the starburst region (primary or secondary) 
are confined to the core but the total possible population 
of electrons/positrons is not equivalent to the confined 
population. This is one of the mechanisms responsible for 
keeping the galaxy on the FIR-radio correlation. With- 
out some loss of protons to a advection in a galactic wind, 
there would be a larger population of electrons/positrons 
resulting in a higher radio flux. 

5.3. Comparison with Other Models 

Our approach to modeling cosmic rays in star- 
bur sj^jrakaxie^l^uflji^ 

els tide Cea del Pozo et al.ll2009at IPaglione fe Abrahams! 
12012ft . However, we place an emphasis in our mod- 
els on varying and constraining a variety of possible 
wind (advection) speeds rather than assuming a sin- 
gle, static wind speed (Ide Cea del Pozo et al.l l2009at 
IPaglione fc Abrahams! 120121 ). Additionally, whereas 
some previous mo dels have had more free parameters 
(|Lacki et al.ll2010ft . our approach reduces the varied pa- 
rameters to a select few. All other parameters are chosen 
to be in agreement with observations. Thus, we compare 
large numbers of models and strive to produce fits to 
the radio and 7-ray spectra that are more intuitive and 
o bservationally informed . 

IPaglione fc Ab rahams (|2012f ) favor higher mean gas 
densities, ranging from 100 to 1000 cm~ 3 , and strong 
magnetic fields, ranging from 200 to 1000 fj,G. Their best 
solution resulted in a density of n = 600 cm -3 and a mag- 
netic field strength of B = 450 pG. In contrast, we favor 
smaller magnetic field strengths (225 to 350 /-tG). How- 
ever, though we test smaller densities (280 to 415 cm -3 ), 
our best fit s have a comparable density of 550 cm -3 (see 
Table 1). Ide Cea del Pozo et al.l (|2009al ) favor models 
with larger supernova rates (0.1-0.3 SN yr -1 ) and use 
a significantly smaller average ISM density (180 cm -3 ). 
They produce models with magnetic fields strengths from 
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120 to 290 ^G. This range is comparable to our own re- 
sults. 

Li ke our model, both Ide Cea del Pozo et all (|2009al) 
and iPaglione fe Abrahams! (|2012l ) incorporate free- free 
absorption into their radio models with spectra that turn 
over at ~700 MHz and ~200 MHz respectively. Neither 
incorporates a halo synchrotron emission component to 
allow for c omparison with low frequency radio data. Ad- 
ditionally, IPaglione fe Abrahams! (|2012D assume a con- 
stant convection loss timescale of 1 Myr which is suffi- 
ciently longer than the energy loss timescale for cosmic 
rays (at their preferred average density of 600 cm -3 ) that 
they find the galaxy to be calorimctric for both protons 
and electrons. 

In addition to creating a model for cosmic rays in star- 
burst galaxies, we address the question of calorimetry in 
starburst ga laxies. Calorimet r y is a lso address by the 
recent paper lAckermann et al.l ( 20121) . A semi-analytical 
formula which depends on the star-formation rate and 
supernova acceleration efficiency of a galaxy is used to 
determine the calorimctric efficiency estimate, while we 
compare a fully calculated energy loss lifetime against 
the advection t imescale for a variety of wind speeds (sec 
§3.1 and §5.1). lAckermann et all ()2012D find that galax- 
ies with SFR ~ 10 Mq yr^ 1 have calorimctric efficiencies 
of 30% to 50%, a similar result to our own findings. 

6. CONCLUSIONS 

While studying the complex situation presented by cos- 
mic ray interactions in the starburst galaxy M82, we 
found it useful to develop a physically realistic but com- 
putationally simple single zone model. In this paper, we 
present a cosmic ray model for M82 that has few free 
parameters and is observationally informed. Our model 
produces statistically reasonable fits to both the radio 
and 7-ray spectra. We find that the inclusion of sec- 
ondary electrons and positrons as well as a galactic wind 
are key to producing better fits. 

In the course of modeling the radio spectrum, we find 
that geometry is vital due to a rather high ionized gas op- 
tical depth at low frequencies. We achieve the turnover of 
the radio spectrum in the starburst core due to free-free 
absorption, but do not reproduce the observed high fre- 
quency free-free emission. Without adding another level 
of detail, it is impossible to achieve both the spectrum 
turnover at low frequencies due to free-free absorption 



and the necessary flux from free-free emission at high 
frequencies. 

While M82 is undoubtedly an electron calorimeter, it 
is not straightforward to make such a distinction for the 
cosmic-ray protons. With our original assumption of 
an acceleration efficiency of 10%, M82 is a 50% proton 
calorimeter, meaning energy losses within the starburst 
zone equal advection losses in the galactic wind. We 
found that it was not possible to achieve acceptable fits 
for an 80% calorimeter when decreasing the acceleration 
efficiency to 4%. However, of 49,000+ models, fits can 
be achieved for an increase of the cosmic ray accelera- 
tion efficiency to 20% (or with an efficiency of 10% and a 
supernova rate of 0.14 yr _1 ) such that M82 is a 35% pro- 
ton calorimeter. Thus, M82 is at best a 50% calorimeter 
and from these results, it is clear that a wind is vital to 
the modeling of cosmic ray interactions in the starburst 
region. 

One of the implications of the FIR-radio correlation is 
that the magnetic and cosmic-ray energy density scale 
with the SFR. From our models, we find that this is 
partially fulfilled in M82. In particular, our x 2 analysis 
puts the magnetic field near (slightly above) equiparti- 
tion, independent of any minimum energy assumption. 
This is also fulfilled in that we find M82 to be an elec- 
tron calorimeter. Thus, our results go part way toward 
explaining the FIR-radio correlation. However, as we 
do not find the galaxy to be a proton calorimeter, the 
cosmic-ray energy density does not completely scale with 
the SFR. As a consequence of not being a perfect proton 
calorimeter, both the 7-ray and radio spectra are dim- 
mer than they would be otherwise, as the loss of protons 
leads to a loss of secondary pions and thus a loss of both 
7-rays and secondary positrons and electrons. Addition- 
ally, the neutrino flux will be a factor of a few less than 
predicted by calorimeter models. 
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APPENDIX 

Energy loss rates: For protons, losses include ionization, Coulomb effects, and pion production. Ionization losses, 
for a neutral medium, are proportional to the density of the medium and weakly dependent on energy such that 



dE 
~dt 



= 1.82 x 10~ 7 n(l + 0.0185 ln{/3) 9[/3 - 0.01]) 



2/3 2 



Ion,p 



10- 



2/3 3 



eV s" 



(1) 



where (3 = v/c and 9 denotes the heavyside theta function (Schlickeiser 200 2]). In an io nized medium, Coulomb losses 
depend on the electron density (n e ) and temperature (T e ) such that (|Schlick ciscr 2002T) 



dE 
~dt 



3.08 x 10" 7 n t 



Coul,p 



j3 3 + 2.34 x 10- 5 (T e /2 x W 6 K) 3/2 



x e 



/3-7.4x 10~ 



2 x IQ 6 K 



1/2' 



eV s" 



(2) 
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The domi nant loss for proto ns above ~4 GcV is pion production. The threshold proton energy for pion production is 



lie dominant loss tor proto ns above ~l uev is pion production. 
.22 GcV fe chlickeisen[2002T) and the energy loss rate is given by 

1.31 x 1CT 7 n 7 1 ' 28 6[ 7 - 1.3] eVs" 1 . (3) 



dE 
~dt 



Pion,p 



For electrons, losses include ionization, bremsstrahlung, inverse Compton, and synchrotron. Ionization losses are 
weakly dependent on energy and are proportional to density such that 

= 36.18 cct t ji(1 + 0.146 Zn(7) + 1.54 n e ) eV s" 1 . (4) 

V M ) Ione 

where <jt is the Thomson cross section, sigmax = 6.65 x 10~ 25 cm 2 . Bremsstrahlung is the dominant energy loss below 
~few GcV. Energy losses due to bremsstrahlung are proportional to energy and dependent on the medium density 
and the energy loss rate is given by 

dE\ Z 2 e g I 192 \ 

where Ze is the charge of the nucleus (|Longairll2011l ). Inverse Compton s cattering goes as t he square of the electron 
energy and depends on the radiation field energy density, U ra d, such that ([Schlickciscr 2 0021) . 



dE\ 4 / E 



= — C Gt 

dt hc,e 3 



U rad eVB- 1 . (6) 



The synchrotron e nergy loss rate also depends on the square of energy and the magnetic field energy density such that 
(lSchlickeiserll2002l) 

' dE \ 4 ( E V B ' 2 ,r 1 
-77) = o C CT T 2 V~ CV S ■ 7 

dt J Synch.e 3 \ m z C J 2 ^0 

Synchrotron and inverse Compton losses are dominant above a few GcV. 
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